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1. Introduction 

The critical current at which a voltage appears along the length of a superconductor is 
one of the most important factors considered in applications of superconductivity pQ. 
The critical current is known to depend upon both local pinning centers in the material 
and the shape of the conductor's cross section. Even in the absence of bulk pinning, 
isolated type-II superconducting thin-film strips subjected to perpendicular magnetic 
fields show magnetic hysteresis due to geometrical edge barriers El HI EJ EJ [7] . Such 
strips have finite critical currents arising from the edge barriers [Sj. It also has been 
shown that the critical currents become larger when slits are fabricated near the edges 
of the strips [9]; the edge-barrier effects are enhanced because the slits increase the 
number of edges that can prevent flux penetration into the inner strips. Although at 
subcritical applied fields and currents the magnetic response of two parallel strips in 
the Meissner state is reversible [10], when the strips are connected at their ends and 
an applied magnetic field exceeds a certain value, magnetic flux penetrates irreversibly 
and the magnetic response becomes hysteretic as a consequence of the edge barriers 
[TT1 [T2l [T3] . A detailed analysis of the effects of edge barriers upon the magnetization 
hysteresis in samples consisting of one, two, and three parallel strips connected at their 
ends has been presented in [14j . 

To calculate the combined effects of geometrical barriers and bulk pinning is more 
difficult, but this has been accomplished for single strips in [U [51 [7J [151 EE E3 US f!9] . 
and a theoretical analysis of the magnetic-field dependence of the critical current of an 
isolated superconducting strip due to both an edge barrier and uniform bulk pinning 
has been presented in [20J. Here we extend the above calculations to account for both 
geometrical edge barriers and bulk pinning in an infinite number of strips using the 
A-array method |21j . by which the magnetic-field and current-density distributions 
for an array of parallel superconducting strips arranged periodically along the x axis 
in the plane y = can be obtained analytically from the solutions for an isolated 
superconducting strip |20j. In particular, we consider the case for which each strip in 
the array carries an equal amount of current in the presence of a perpendicular applied 
magnetic field. We then calculate the critical current I c of each strip accounting for 
both a geometrical edge barrier, which impedes the entry of vortices into the strip, and 
uniform bulk pinning, which impedes the motion of vortices across the strip. 

Our calculation is relevant to a number of recent studies of the ac properties of 
striated coated-conductor tapes, i.e., tapes that have been subdivided into parallel thin 
strips (filaments) separated by narrow gaps. Ideally, to minimize ac losses in power 
engineering applications using multifilamentary conductors, the individual filaments 
should not only have small cross sections but also be twisted and transposed, such 
that the filaments are inductively equivalent, are decoupled from each other, and 
thus carry equal current [22J. In practice, it is not possible to satisfy all these 
conditions, and compromises are generally accepted for practical reasons. Carr and 
Oberly [23J have suggested that that the ac losses in tapes several millimeters wide 
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Figure 1. Long thin superconducting strip of thickness d and width 2w. 



could be reduced by first subdividing the tapes by striations and then twisting the 
tapes. Analytical and experimental results pursuing this idea have been reported in 
[2H EHl [261 EZl EHl [301 EH E2l [331 EH EHl [361 EZJ- Ashworth and Grilli [38] have 
recently proposed a variation of this approach, in which the tape is also subdivided into 
narrow parallel filaments, but instead of twisting the tape, the filaments are interrupted 
periodically by transverse cross-cuts bridged with normal metal. Magnetic flux can 
enter at the cross-cuts, thereby decoupling the filaments and allowing equal currents to 
flow in each filament. The additional ohmic losses in the normal bridges are more than 
compensated by a huge reduction in the far more important ac coupling losses. 

In the present paper we examine the dc properties of an infinite array of parallel 
superconducting strips, which can be regarded as an approximation to a striated coated- 
conductor tape of finite width. Our paper focuses on the possibility that a geometrical 
edge barrier can increase the critical current of each strip above that due to bulk pinning 
alone, but it also should be possible to use our results in calculating the hysteretic ac 
losses of an array of parallel strips as in [211 ESI El IH] • 

Our paper is organized as follows. In Sec. [2] we review our complex-field approach 
and the A-array method. In Sec. [3] we apply the complex-field approach and the A- 
array method to obtain the complex field and the critical currents in a periodic array of 
parallel superconducting strips subject only to bulk pinning. In Sec. @]we calculate the 
magnetic-field dependence of the critical current of one of the strips when the infinite 
array is subject only to geometrical edge barriers. In Sec. [5] we calculate the critical 
current of one of the strips when both geometrical edge barriers and bulk pinning are 
present. We discuss the results and present our conclusions in Sec. El 

2. Complex fields and A arrays 

We begin by reviewing the properties of a long superconducting strip of thickness d 
centered on the z axis in the region \x\ < w and \y\ < d/2, where d <C w [see figure 1]. 
(We attach tildes to all quantities related to the single strip.) Since we are not concerned 
with details of how the current density varies across the film thickness, we consider only 
the sheet-current density K(x) = J_ d , 2 jdy = zK z (x). If the film thickness d is less 
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than the London penetration depth A, we assume that the two-dimensional screening 
length (or Pearl length) [12] A = 2X 2 /d is much smaller than w. We consider the 
general case when the strip carries a current in the z direction and is subjected to a 
perpendicular applied field H a in the y direction. Outside the strip, the net magnetic 
field H(x, y) = xH x (x, y) + yH y {x, y) obeys V x H = and V • H — 0. 

For a two-dimensional problem such as this, it is convenient to introduce a complex 
magnetic field H(C) = H y (x, y) + iH x (x, y), which is an analytic function of the complex 
variable ( = x + iy outside the strip. The real and imaginary parts of H(Q obey the 
Cauchy relations, dH y /dx = dH x /dy and dH y /dy = —dH x /dx, which guarantee that 
H obeys V x H = and V • H = 0, respectively. The Biot-Savart law can be expressed 
as 

H(()^.4f^. (1) 

^ J-w (-U 

Using the property that l/(x ±ie) = {P/x) =F i5(x), where e is a positive infinitessimal 
and P denotes the principal value, we obtain 

^»)=i + fr^ ( 2 , 

2?r x-u 

and 

H x (x, ±e) = =F K z (x)/2, \x\ < w, 

= 0, otherwise. (3) 

The net current carried in the z direction is 

I z = / K z (x)dx. (4) 



Note from ([T]) that an expansion of H(() in powers of l/( yields [H] H(Q = 
H a + V2vrC + 0(l/C 2 ). 

Solutions for the magnetic-field and current- density distributions for a long thin 
strip are known for many different physical situations. It has been shown by Mawatari 
in [21] how known solutions H{() and K z (x) for a single isolated strip can be used to 
generate the corresponding solutions H(() and K z (x) for an X array, i.e., an infinite 
periodic array of identical coplanar strips of width 2w and periodicity L in the plane 
y = 0, as sketched in figure 2. The complex field H(() = H y (x,y) + iH x (x,y), which 
is an analytic function of the complex variable ( = x + iy outside the strips, describes 
the components of the net magnetic field H(x,y) = xH x (x,y) + yH y (x,y) produced 
by a sheet-current density K(x) = zK z (x) and a magnetic field H a applied in the y 
direction. The periodicity requires that the magnetic field and sheet-current density 
obey H(x + L,y) = H(x, y) and K{x + L) = K(x). 

Corresponding solutions for the periodic complex field H(() can be obtained from 
H(() using the conformal mapping [21] 

C = (L/7r)tan(7rC/L), (5) 
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Figure 2. X array of long thin superconducting strips of thickness d, width 2w, 
periodicity L, and separation A = L — 2w. 



or its inverse 

C = (L/tc) arctan(7rC/L). (6) 

These equations also describe the relations between x and x, w and w, or u and u, 
which will appear in later expressions. Then H(() = H(C), K z (x) = K z (x), and the 
Biot-Savart law yields [21] 

i r .. , . ^(c-u)- 



H(() = H a + 



2L 



KJu) cot 



du. 



where the magnetic field applied in the y direction in the £ plane is [21 

1 P K z {u)udu 



H a = H a 



2tt J_a (L/tt) 2 + v? 

and the current carried in the z direction by one of the strips shown in figure [2] is 

K z (x)dx. 



(7) 



(8) 



(9) 



Note by comparing (j4j) and (Q that I z is in general not the same as I z . The relations 
dHJ) and (jHJ) can be obtained from (JTj) by applying the requirement that 

H{±ioa) = H(±iL/n) =H a T U,/2L. (10) 

Along the x axis, we have the properties H(x±ie) = H(x±ie), H x (x±ie) = =pK z (x)/2, 
and 

- / ^(n)tan( T 



H y (±L/2, 0) = H a = H a + I K z (u) tan(^f)dw. 



'11' 



In the following sections, we shall obtain H{Q using i?(C) = H(Q, K z (x) = K z (x) 
(jSJ) and pop , and various trigonometric identities. 



3. Critical current due only to bulk pinning 



When the critical current of a single superconducting strip is due solely to bulk pinning, 
characterized by the bulk-pinning critical current density J p (here assumed to be field- 
independent, an assumption to be justified in Sec. ED, the corresponding critical sheet- 
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current density is K p = J p d, and K z = K p for \x\ < w, such that the critical current is 
I p = 2K p w and the complex field is 



g(C) = g. + ^ln r[^ +m i^ -, (13) 



The corresponding complex field surrounding the X array shown in figure [2] is 

-sm [vr(C-w)/L]. 

where H a = H a = H y (±L/2, 0). Since K z = K z = K p , we immediately find for one 
of the strips in the X array that its critical current is I p = 2K p w, its average critical 
sheet-current density is K p , and its average critical current density is J p . Note that when 
w < L/2, the self-field contribution to H y (x,0) is positive and diverges logarithmically 
when x — > w; it is negative and has a similar divergence when x — > —w. In the limit 
that w — > L/2, however, we recover the complex potential generated by an infinite film 
carrying a sheet current K z = K p in a perpendicular applied field H a , 

H(() = H aT iK p /2, (14) 

where the upper (lower) sign holds when y > (y < 0). 

4. Critical current due only to edge barriers 

Consider a single bulk-pinning-free superconducting strip in which the critical current is 
determined by a geometrical edge barrier (or surface barrier at the edge). As discussed 
in [5] , the current and field distributions at the critical current have two possible forms, 
depending upon the strength of the applied field H a . For small values of H a , the 
distributions are simply the Meissner response to the applied field and current, but for 
larger values of H a , there is a domelike magnetic field distribution inside the strip. 

The complex magnetic field describing the Meissner response to a magnetic field 
H a applied in the y direction and to a net current / applied in the z direction is [5] 

H(0 = ¥ + I ' 2 \ (15) 
and the sheet-current density in the strip (|x| < w) is 

~ 2H a X + I/7T 

K z {x) = — (16) 

[w — X ) 1 

If the strip contains a domelike magnetic-flux distribution in the region a < x < b, 
where the sheet-current density is zero, the complex magnetic field is [5] 

^V^K^Hfz^, (17) 

L Q z — w z J 
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KJx) 



, (x — a)(x — b) 
2H a \l K _J y _„ \ b<x<w, 



w 2 — X 2 



(a — x)(b — x) 
w 2 — x 2 



— w < x < a, 



(19) 



= 0, otherwise, 
and the net current in the z direction is [5] 
J = -7r(a + b)H a . 

The corresponding quantities for the X array are, for the Meissner response to 
an applied field H a in the y direction and an applied current / in the z direction, the 
complex field 

(20) 



tf a sin($) + (//2L)cos( L 



[ si n(^±^) s in(^ r ^)]V2 

and the sheet-current density in the strip (\x\ < w 

2H a sm(f) + (I/L) cos 



KJx) 



nx ' 
L ' 



sm (ZL(^h£)) sin (^£l; 



(21) 



where H a = H a / cos(nw/L) = H y (±L/2,0) and 1 = 1/ cos(ttw/L). 

If all the strips in the X array contain domelike magnetic-flux distributions identical 
to the one in the region a < x < b, the complex field is 



H(0 



cos 



- 7r(a+6) > 
• 2L - 



sinf^F^W 7 ^^ 



L) sin (ZLikz£2) , 1/2 



sin (ZL(^) s i n (l(C^); 



(22) 



the sheet-current density in the central strip (\x\ < w) is 



KJx) 



2H a 



cos 



• ir(a+b) ■ 
■ 2L ■ 



2H n 



sin 



■ ir(x-a) 



) sin 



ir(x — b) ' 
I - 



\ ■ I ir(w+x) \ ■ i 

\ sml^J sm( 



tt(w— x) ^ ' 

b < X < W , 



COS 



ir(a+b) 
2L 



A 



sin 



7r(a— x) 



) sin( 



7r(6— x) 1 



sin -H- — '- sin -H — '- 



w < x < a, 



0, otherwise, 



where 



H. 



aV /cos(f)cos(f) 
cos(^)cos(^)" 



tf,(±L/2,0), 



and the net current in the z direction in the strip \x\ < w is 
I = -2LH a tan[7r(a + b) /2L\ . 



(23) 



(24) 



(25) 
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The divergences in the above expressions for H and K z at x = ±u> are artifacts of 
ignoring the finite thickness d of the film. It is well known that these divergences are 
cut off at a length scale A c , the larger of A = 2X 2 /d and d/2 [i.e., A c = max(A, d/2)]. 
To determine the critical current of one of the strips of the X array, we use the 
approximations applied in [20J and [13]. Accordingly, we estimate the the sheet-current 
density at an edge by evaluating the diverging inverse square root in the expression for 
K z at a distance A c from that edge; for example, for the edge at x = w, in (|2T]) and fl23l) 
we replace x by w in the numerators and by x c = w — A c in the denominators and use 
the fact that the cut-off length scale obeys A c C w. To account for the edge barrier, 
we assume that vortices nucleate and enter the superconductor when the magnitude of 
K z at either edge of the strip reaches the critical value K s = j s d at which the barrier 
is overcome. For an ideal edge, j s is equal to the Ginzburg-Landau depairing current 
density jcL [El EH], but for an extremely defected edge, j s may become very small. 
Applying this procedure to (I2T1) for H a = 0, we obtain the following approximation to 
the zero-field surface-barrier critical current for any one of the strips in the X array, 

I s (0) = I s0 = K s (2vr A C L tan 9) V 2 , (26) 

where 9 = ttw/L. Alternatively, by estimating the local magnetic field at the edge of 
the strip [51 S3], one obtains ( |26i) but with K s rs 2H s , where H s is the critical field at 
which the barrier is overcome. We expect that H s is at least as large as the lower critical 
field H c i, but under favorable circumstances it may approach the bulk thermodynamic 
critical field H c . 

As a function of the applied field H a , the surface-barrier critical current is found 
from ( T2T|) to obey 

I,(H a )/I a (0) = l-h (27) 
for small h, where 

h = H a /(I s0 /2Ltan9). (28) 

When L — ► oo, ff26l) and (T27j) reduce to corresponding results found in [5] and (20] in 
low fields for isolated bulk-pinning-free strips. 

In the linear region of I s (H a ) given in (1271) . the flux flow producing the voltage at 
a current just above I s can be described as the nucleation of vortices at x = w as they 
overcome the geometrical edge barrier, followed by rapid motion across the strip and 
annihilation at x — —w. This can occur only when K z (x) > for all \x\ < w. For 
increasing H a , however, as can be seen from ( l2~Tj) . K z (—w) becomes zero at the critical 
current I s {H a ) when H a = I s (H a )/2Lt&n9. Combining this condition with (l2"Tj) . we 
find that this occurs when H a = H^o, where 

H d0 = I s0 /4Ltan9, (29) 

and the subscripts (dO) denote the onset of a dome for zero bulk pinning. Note that 
we also may write h = H a /2Hdo- The linear portion of I s (H a ) vs H a given in ( 1271) ends 
when H a = H d0 or h = h d0 = 1/2. 
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For applied fields H a > H dQ (or h > h d0 = 1/2), each strip contains a domelike 
magnetic field distribution at the critical current, and I s (H a ) is no longer a linear 
function of H a . The flux flow producing the voltage at a current just above I s can 
be described as the nucleation of vortices at x = w, as they overcome the geometrical 
edge barrier, followed by rapid motion from x = w to x = b, the right boundary of 
the magnetic field dome. The vortices inside the dome move collectively very slowly, 
and vortices at x = a « — w get pushed out of the dome and annihilate at x — —w. 
We therefore calculate the critical current by setting a = —w in ( 123|) . evaluating K z at 
x = w using the above approximation method, and setting it equal to K s . Using (1261) 
and ( 128]) we then find that the value of h at the critical current is b c , where 

tan[vr(> - b c )/2L] = tan(vrw/L)/4/i 2 . (30) 

(Note that b c = —w when h = 1/2.) Substituting this into (125]) . we obtain I s (H a ) for 
H a > Hdo (or h > hdo = 1/2). The critical current of one of the strips of the X array is 
thus given by 



, h, 0<h< h d0 

'sO ^ 



— 7T) hdo — ^ < h < h max , (31) 
Ah 2 

when there is no bulk pinning and the critical current is determined only by a geometrical 
edge barrier (surface barrier). We expect our theory to be accurate only at sufficiently 
low applied fields that h < h max or H a < H max , where h max is the reduced field at which 
b c = w - A c . From ([30]) we obtain h max = [(L/2ttA c ) tan^] 1 / 2 > (u>/2A c ) 1/2 > 1 and 

±J -max ~ 

5. Critical current due to both edge barriers and bulk pinning 

In this section we calculate the complex fields and critical current J c of one of the 
strips in the X array when the critical current is due not only to a geometrical edge 
barrier but also bulk pinning, characterized by a bulk-pinning field-independent critical 
sheet-current density K p . An important parameter is the ratic[jj 

ip sine 

p ~i s0 e ' {d2) 

where I p = 2K p w [see Sec. [3] is the critical current due to bulk pinning in the absence 
of an edge barrier, 7 s0 [see Sec. H] is the critical current due to the edge barrier in the 
absence of bulk pinning when H a = 0, and 6 = irw/L. 

5.1. Region I: One vortex-free zone with no domes 

We first consider the case of relatively weak bulk pinning when p < 2/n. In low fields 
(0 < H a < H d or < h < ha, region I of figure 3), vortices nucleate at the right edge 

|| In the limit L — > oo, p = I p /I s o as in [2U] . 
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Figure 3. Diagram showing the behavior at the critical current vs reduced field h 
and bulk pinning parameter p. In region I, the strip is vortex- free [H y (x,0) — 0]. In 
region II, there are two zones: a vortex dome (— w < x < b) and a vortex-free zone 
(b < x < w). In region III, there are three zones: a vortex dome (a < x < b) between 
two vortex-free zones (—w < x < a and b < x < w). In region IV, there are four 
zones: an antivortex dome (a < x < c) where H y (x, 0) < 0, a vortex dome (c < x < b) 
where H y (x, 0) > 0, and two vortex-free zones (— w < x < a and b < x < w). In each 
vortex-free zone K z (x) > K p , and under each vortex dome K z = K c . The curve hd(p) 
separates regions I and III, the curves hi (p, 6) separate regions III and IV, and the 
curves /i2(p, 9) separate regions II and III. These curves are shown for 9 = ixwjL = 
(solid), 6 = 0.5tt/2 (dot-dashed), and 9 = 0.99?r/2 (dashed). 



of the strip (x = w) when I z slightly exceeds I s (H a ). When K z (x) in (12TT) is greater 
than K p for all x in the strip (|a;| < w), nucleating vortices are driven entirely across 
the strip to the opposite side at x = —w, where they annihilate. The critical current is 
then I c (H a ,p,6) = I s (H a ), and the normalized critical current is 



the same result as in (|27|) . Thus if p is small, the critical current I c (H a ,p,8) for small 
H a is still completely dominated by the edge barrier and is independent of the strength 
of the bulk pinning. From (1211) we see that K z (x) vs x has a minimum, K z (x m i n ), at 
x = x min , where 



The minimum deepens [i.e., K z (x iain ) decreases] with increasing H a , and K z (x min ) = K p 
when / = I c (H a ,p, 6) and H a = Ha or h = ha, such that 




(33) 




(34) 



tan(7rx min /L) 



[h d /(l-h d )}ttm6, 



(35) 
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where 



h„ 



1 



7T P^2^ 



-[l-( — 

2 L v 2 



(36) 



I s0 /2Ltane 

As a result, for small p < 2/tt, the linear behavior of I c (H a ) vs H a , given by i c (h,p,9) 
vs h in (l33l . holds only for < h < hd. 

When p < 2/n and h > hd, or when p > 2/n and h has any value, i.e., for h and 
p outside region I of figure 3, domelike magnetic- field distributions occur at the critical 
current I c (H a ). Using the X-array method [21], we can obtain the complex field H(() 
and associated sheet-current density K z (x) in one of the strips of the X array from the 
corresponding solution for an isolated strip [20j . For the X array, we find H(Q = H((), 
where 



H(C) = P(a,b,C)[H a + ^Q(a,b,C)], 

In 



where 



P(a,b,C) 



(C-fi)(C-&) 

(C 2 - w 2 ) 



1/2 



(37) 



(38) 



Q(a,b,() 



and 



\fw 2 — u 2 



(C- 

2(w 



-du 



u 



a)(b 



u 



(w — a){w + b) 



n 



b — a 



n 



■w + b 
w){b 



(C-5) \((-a)(w + b) 



a) 



2(w - b) 



n 



w 



(C 



a)(w + 6) 
(C-«)(6 



b — a 



w 



n 



(C-6) HC-6)(t3-o) 



? 2 



2w(6 - a) 



(39) 



(40) 



(41) 



(42) 



(w — a)(u; + 6) 

The integral in (1391) is expressed in terms of complete elliptic integrals of the third kind 
[4"6| Wf[ HH1 H9j [50] II(n, /c), where n is called either the characteristic or the parameter, 
and k is called the modulus. The mappings of ^ and ^ define the relations between 
( and (, a and a, b and b, or w and w, and also guarantee that H(() is periodic in the 
x direction with period L. The applied magnetic field H a in the y direction and the 
current I z carried in the z direction by one of the strips shown in figure [21 obtained from 
(HO]), obey 



Re{P(5, 6, iL/n) [H a + ^Q(a, b, iL/ti)}} 

Z7T 



(43) 
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and 

I z = -2LIm{P(a,b,iL/n)[H a + ^Q(a,b,iL/7r))}, (44) 

where 

P(a, b, lL /tt) = COS9 e ™(a+b)/2L (45) 

\J cos(7ra/L) cos(irb/L) 

and 9 = irw / L. 

The quantities of primary interest to us are H y (x,0) = ReH(x) and K z (x) = 
-2lmH(x + ie), 

H y (x,0) = ReH(x) (46) 

= P {a,b,x)[H a + -^ReQ{a,b,x)}, (47) 

a < x < b or \x\ > w, 
= 0, —w <x<aoib<x<w. (48) 

K z (x) = K z (x) = -21mH{x + it) (49) 
= -2P (a,b 1 x)[H a + ^Q(a,b,x)}, (50) 

Z7T 

— w < x < a, 

= K p , a < x < b, (51) 

= +2P (a,b,x)[H a +^Q(a,b,x)], (52) 

Ztx 



b < x < w, 



where 



\x 


— a\ 


\x — 


b\ 




x 2 - 


- w 2 \ 





P (a, b,x) = J L^L-^. (53) 



For the calculation of H y (x, 0) for a < x < b (i.e., a < x < b), taking the real part of Q 
in (T4T|) corresponds to taking the principal value of the integral of (1391) . 



5.2. Region IV: Two vortex-free zones, one vortex dome, and one antivortex dome 

For values of p > 2/n [see (1321) ]. and small values of h, i.e., for h and p in region 
IV of figure [HI the vortex distribution in each of the strips at the critical current can 
be described as a double dome, consisting of a vortex dome adjacent to an antivortex 
dome. Just above the critical current, vortices nucleate at x = w, where K z (x c ) = K s , 
move rapidly to the left through an otherwise vortex- free region (6 < x < w), and then 
move slowly to the left through a vortex-filled region (the vortex dome), c < x < b. 
Antivortices nucleate at x = —w, where K z (—x c ) = K s , move rapidly to the right 
through an otherwise vortex- free region (— w < x < a), and then move slowly to the 
right through an antivortex-filled region (the antivortex dome), a < x < c. Vortices and 
antivortices annihilate at x — c, where the two domes meet and H y (c, 0) = 0. 
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Four equations must be solved simultaneously for a, b, H a} and the critical current 
I c (H a ) for known values of h, p, and 6 = ttw/L in region IV of figure [3) One condition 
is that K z (w) = K s , which we evaluate by replacing x in the denominator of ( |53i) by x c , 
where x c = w — A c and x c = (L/tt) tan(7nr c /L) rs w — A c sec 2 9, and by replacing x by 
w in Q(a, b, x) (HD1 and in the numerator of ( 1531) . This yields from (132]) and (152]) 

yJ(l-a')(l-V)hcQs6 + (1 + a')7^n(^i, g> = 1, (54) 

where we use the normalized quantities a' = a/w and 6' = b/w, as well as the definition 

h = H a /(I s0 /2Lt&n9). (55) 

A second condition is that K z (—w) = K s , which we evaluate in a similar manner with 
the help of (l3"2l), (jiljl. and CT. The result is 

- \/([l + a'){l + b')hcos9 + (1 - 6')\ /i±^n(^^-,g)p = 1. (56) 
v V 1 — cl 1 — a 

Equations fl54l) and fl56|) have nearly the same form as (9) and (10) in [20J and reduce 

exactly to those equations in the limit as 6 = ttw/L — > or L — > oo. The third and 

fourth equations needed are obtained from (143]) and (H4l . which we write in dimensionless 

form at I z = I c as 

h=P 1 h + (p/2)(AQi - P2Q2) sec 0, (57) 
i c = - cot 6[P 2 h + (p / 2) (P 1 Q 2 + P 2 Q 1 ) sec 6], (58) 

where the real quantities Pi,P 2 ,Qi and Q2 are defined via P(a,b,iL/n) — Pi + «A 
and Q(a,b,iL/ir) — Qi + iQ 2 , h, h, and p are defined via (|28l) . ( J55l) . and (|32l) . and 
i c (h,p,9) = I c /I s o- F° r given values of the strip width 2u>, periodicity length L, 
dimensionless applied field h, and dimensionless bulk pinning strength p, numerical 
solutions of ( 1541) . ( |56|) . ( |57|) . and ( 1581) yield the strip's dimensionless critical current 
i c (h,p, 9) as well as the three other unknowns, a (or a), 6 (or b), and /i. Shown in figure 
H] (a) and (b) are sample plots of H y (x, 0) and K z (x) for /i and p in region IV. 



5.3. h = 



Simplifications occur when calculating the self- field critical current when p > 2/n and 
H a = {h = 0), for which a = — b (a = —6) and H y (x, 0) is antisymmetric about the 
point x = c = 0. Equations (138]) and (139]) then can be evaluated as 

b 2 - C 2 \ V2 



w 2 — ( 2 



Q(-6,6,C) = 2i 



K 



, (^ 2 -c 2 ) n /& 2 1 

C 2 V C ~2'^yj 



(59) 

(60) 



where K(fc) is the complete elliptic integral of the first kind of modulus k. In 
and (jB]), P(—b,b,iL/ir) is pure real and Q(—b,b,iL/7i) is pure imaginary, such that 
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Figure 4. Plots of (a) H y (x,0) (in units of I sO /2Lta,n0) and (b) if z (x) (in units of 
K p ) vs a; (in units of w) for h — 0.2, p = 0.85, and 9 = ttw/L = 0.57r/2 (region IV), 
for which i c = I c /I s q = 1-11 a-nd / c // p = 1.18, showing three strips (—5 < x/w < —3, 
— 1 < x/w < 1, and 3 < x/w < 5) and the gaps around them (—7 < x/w < —5, 
-3 < x/w < —1, 1 < x/w < 3, and 5 < x/w < 7). 



from (|57j) we obtain h = when /i = 0. Equations (|54l) and ( 1561) . resulting from the 
requirement that K z (w) = K z (—w) = K s at the critical current, reduce to a single 
equation, 



P K{b')Vl-b' 2 = 1, (61) 

and the expression corresponding to fl58l) for the reduced critical current i c becomes, 
with the help of (ED), 



' //2) -nf - & 1 . 9n ,b')/K(b% (62) 



(l + 6 /2 tan 2 ^) Vcos 2 + 6 /2 sin^ 

where 9 = ttw/L. When b' = 0, we find that p = 2/ir, and i c = 1 for any value of 9. 
When = 0, (El reduces to 



z c = . V-^ - (63) 
K(6')vl -&' 2 

Plots of H y (x,0) and vs x for h = are similar to figure H] (a) and (b), except 

that H y (x,0) is an antisymmetric function of x, centered at the origin, and K z (x) is a 
symmetric function of x. 
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5-4- hi, the boundary between regions III and IV 

For increasing values of the applied field H a (or the reduced field h), the vortex dome 
expands and the antivortex dome shrinks. For reduced fields h in the range < h < hi, 
we have H y (x,0) < for x slightly larger than a. In other words, the coefficient of 
Po{a, b, a + e) in (T4T1) is negative. However, when h = hi, this coefficient becomes zero. 
For h > hi, this coefficient is positive, indicating that at the critical current, only a 
vortex dome is present (region III in figure [3]). The value of hi is determined chiefly by 
the condition that the coefficient of Po{a, b, a + e) in (H7|) is zero, which yields with the 
help of flUD, (17.7.7) of g7], and (414.01) of [50], 



h + psec6 [(i - a')(i + mm 



where 



(b ! -a')^(l-a')(l + V) 

-(l + ~a')(l-V)K(q") 

_ { y - a')(l - b')U( b ^4,q')} = 0, (64) 
1 — a' J 

,9 2(1' -a') , s 

q' 2 = i 65 

(l-5')(l+&0 

The value of hi for given values of p > 2 /it and 9 = ttw/L is obtained as the value of h 
that satisfies fl64l) . fl54|) . and fl56l) . The solutions of these equations also yield the values 
of a' and V at the critical current when h = hi. 



5.5. Region III: Two vortex-free zones and one vortex dome 

In region III of the diagram of p vs h shown in figure [3], at the critical current, the field 
and current distributions within each strip divide into three zones. For b < x < w, 
we have H y (x,0) = and K z (x) > K c ; for a < x < b, we have H y (x,0) > and 
K z (x) = K c ; and for — w < x < a, we have we have H y (x,0) = and K z (x) > K c . 
Just above the critical current I c , the voltage along the length of the strip is produced 
by (a) small numbers of vortices nucleated at x = w as they overcome the geometrical 
edge barrier, followed by their rapid motion from x = w to x = b, (b) large numbers 
of vortices moving slowly through the vortex dome from x = b to x = a, and (c) small 
numbers of vortices leaking out of the vortex dome and moving rapidly from x = a to 
x = —w. 

The four equations determining the values of a, b, h, and the reduced critical current 
i c for given values of h, p and 9 in region III are (|5T|) . ([SHD, and (1641) . since the 

latter equation also can be shown to give the condition that dK z (x)/dx = at x = a. 
Shown in figure (a) and (b) are sample plots of H y (x, 0) and K z (x) for h and p in 
region III. 
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Figure 5. Plots of (a) H y (x,0) (in units of J s o/2Ztan0) and (b) K z {x) (in units of 
if p ) vs cc (in units of w) for ft, = 0.2, p — 0.6, and = irw/L = 0.57r/2 (region III), 
for which i c = I c /I s q = 0.87 and I c /I p = 1.31, showing three strips (—5 < x/w < —3, 
— 1 < x/w < 1, and 3 < x/w < 5) and the gaps around them (—7 < x/w < —5, 
-3 < x/w < —1, 1 < x/w < 3, and 5 < x/w < 7). 

5.6. ti2, the boundary between regions III and II 

For increasing values of h, the left boundary of the vortex dome moves ever closer to 
the left edge of the strip, i.e., a — > w or a — > —w, and numerical difficulties arise. 
Although other criteria could be used, in this paper we state that for practical purposes 
the vortex-free zone on the left side shrinks to negligible width when a = — 0.9999w. 
For given values of p and 9, this occurs at a reduced field h = h^, which can be obtained 
by solving ( |54l) . ( |64l) . and ( |57|) for h, h, and b after replacing a by — 0.9999w. 

5. 7. Region II: One vortex-free zone and one vortex dome 

For values of h > h 2 in figure [3j all the quantities at the critical current can in principal 
be calculated using the same equations as for Region III. However, these quantities can 
be calculated with fewer numerical difficulties and with high accuracy by setting a = —w 
in the above equations, which results in a number of simplifications. The complex field 
still obeys H(Q = H(Q, but now, to good approximation, 

H(C) = P(-w,b,()[H a +^Q(-w,b,C)], (66) 

Z7T 
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where 



and 



C - b \ V2 

P(-w,b,0=[^-) (67) 



Q(-w,b,()=[ —^^ =d u (68) 
J-w [C, — u)y b — u 

= 2 sinh- fi±i) "i 2(^3) 1 ^h-f(g±M=^l 1/2 (.69) 

This equation is equivalent to (14b) of [20], which was misprinted without the prefactor 
before the second sinh -1 term. 

H y (x, 0) = KeH(x) and K z (x) = —2lmH(x + ie) in region II are given by 

H y (x,0) = ReH(x) (70) 

= P (-w, b, x) [H a + ^ReQ(-w, b, x)], (71) 

— w < x < b or \x\ > w, 
= 0, b < x < w. (72) 

K z (x) = K z (x) = -2lmH(x + ie) (73) 
= K p , -w < x < b, (74) 

= + 2P (-w, b, x)[H a + ^Q(-w, b,x)), b<x<w, (75) 

Z7T 

where 



Po(-ir.b.:r) = i/J?—^. (76) 

|x — w\ 

Three equations must be solved simultaneously for b (or b), h, and i c for given 
values of h, p, and 9 in region II of figure [31 The condition that K z (w) = K s , evaluated 
as in deriving (1541 . becomes 



2(1 -b>) ( /icostf + psinh- 1 \l^-j t J = X > ( 77 ) 

and the other two equations needed are (!57|) and (158]) . in which the real quantities 
P\,P2,Qi and Q2 are defined via P(—w,b,iL/n) = P\ + iP 2 and Q(—w,b,iL/n) = 
Qi + Shown in figure [6] (a) and (b) are sample plots of H y (x, 0) and K z (x) for 
and p in region II. 



5.5. Critical current 

Using the above approach, one can calculate the normalized critical current i c = I c /I s o 
( T2S1 for arbitrary values of the reduced field h (1281) . pinning parameter p (1321 . and 
9 = irw/L, where the strip width is 2w and the periodicity length is L. Shown in figure 
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Figure 6. Plots of (a) H y (x,0) (in units of I sO /2Lta,n0) and (b) K z (x) (in units of 
K p ) vs x (in units of w) for h — 1.0, p — 0.15, and 9 = ttw/L = 0.57r/2 (region II), 
for which i c = I c /I s q = 0.35 and I c /I p = 2.08, showing three strips (—5 < x/w < —3, 
— 1 < x/w < 1, and 3 < x/w < 5) and the gaps around them (—7 < x/w < —5, 
-3 < x/w < —1, 1 < x/w < 3, and 5 < x/w < 7). 



Hare sample results for 9 = 0.57r/2, when the widths of the strips and the gaps are 
equal. Because of the effects of the geometrical barrier, at sufficiently low applied fields 
and small values of p, the overall critical current I c can be considerably enhanced above 
the critical current I p = 2wK p due to bulk pinning alone. This occurs because the 
sheet-current density K z (x) in the vortex-free zones near the sample edges can carry a 
supercurrent with a density well above K p . However, in high applied fields the cross 
section of the strip becomes nearly completely filled with vortices, and the sheet-current 
density in the vortex- filled zones cannot exceed K p . 



6. Conclusions 



In this paper we have considered an infinite array (X array) of coplanar equally spaced 
superconducting strips and have calculated the critical current of each strip in the 
presence of a perpendicular applied field, taking into account not only geometrical 
barriers at the edges of the strips but also bulk pinning characterized by a field- 
independent critical sheet-current density K p . We have carried out these calculations 
using the X-array method of [21], which enabled us to use solutions from [20] for an 
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0.5 1 1.5 2 

h 

Figure 7. i c (h,p,8) (critical current normalized to I s q) vs reduced field h for 
9 = 0.57r/2 and various values of the bulk pinning parameter p. The solid curves 
exhibit i c vs h, the dashed curve shows i c at h = hi(p, 9), and the dotted curve shows 
i c at h = h 2 (p,9). For p < 2/ir, i c decreases linearly with h ([55)1 up to hd and then 
decreases more slowly in regions III and II. One curve shows i c for the special case of 
p = 2/tt. For p > 2/tt, i c increases by a few percent in the double-dome region IV 
and then decreases more gradually in regions III and II. In all cases, i c asymptotically 
approaches p9 j sin 9 for large h (short dotted lines along the right side of the figure). 

isolated strip. 

Geometrical barriers at the edges of the strips enhance the critical current above 
what it would be if it were due to bulk pinning alone. At low fields, these barriers delay 
the entrance of vortices and permit the strip edges to remain vortex-free and carry a 
sheet-current density well in excess of K p . Any vortices and antivortices entering the 
strips arrange themselves into domelike distributions, and under these domes the sheet- 
current density remains at K p . As the perpendicular applied field increases, more and 
more vortices are forced into the strip, causing the vortex domes to expand and the 
high-current-density vortex-free zones to shrink, such that the overall critical current 
approaches that due to bulk pinning alone, I p = 2wK p . 

The effects predicted in this paper for parallel arrays of narrow superconducting 
strips would be most easily observed experimentally using magneto-optical or scanning 
techniques in materials with low bulk pinning, such as Bi-2212 (Bi 2 Sr 2 CaCu20 8+ 5), 
where geometrical-barrier effects were first observed [31 H], or in low-pinning type-I 
superconductors, such as Pb, where magnetic flux domes consisting of the intermediate 
state have been observed [51]. However, it is important also to consider the possibility 
of observing these effects in materials with strong bulk pinning. 

As a practical application of the above theory to coated conductors made 
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of superconducting YBCO (YBa2Cu 3 07_,5), let us calculate the critical-current 
enhancement due to geometrical barriers for an isolated strip and then estimate how this 
enhancement is affected by striations. We first consider a long flat strip of thickness d = 
1 ixm. and width 2w = 5 mm. We assume that its self-field bulk-pinning critical current 
density is J p = 1 MA/cm 2 (10 10 A/m 2 ), such that the bulk-pinnng critical sheet-current 
density is K p = J p d = 10 4 A/m and the bulk-pinning critical current is I p = 2wK p = 50 
A. We use (1261) with K s = 2H S , A c = d/2, and L —>■ oo to estimate the zero-field 
surface-barrier critical current I s o- We use the experimental results of [52] as discussed 
in [53], from which we estimate that H c \ = 180 Oe = 1.4 xlO 4 A/m and H c = 3.2 kOe 
= 2.6 xlO 5 A/m at 77 K. We use the conservative estimate that H s = H c \, which yields 
J s o = 27T H c i\fdw = 4.5 A. The corresponding scaling field 

H sca i e = I s0 /2L tan 9, (78) 

which appears in the denominator of fT2"8"l) and f[55|) . becomes H sca i e = H c i^Jd/w = 3.6 
Oe = 286 A/m. The corresponding value of the pinning parameter is p = 11, for which 
we obtain from (16T1) and ( 1631) that the self-field critical current I c at h = H a / H sca i e = 
is very nearly equal to I p . Moreover, I c remains very nearly equal to to I p as the 
applied field increases. Our solutions show that the vortex-free regions are very close 
to the edges at the critical current (i.e., b and —a are nearly equal to w), such that the 
geometrical-barrier enhancement of I c is negligible. 

Next we suppose that the 5 mm strip is subdivided into 50 parallel strips, each 
of width 2w = 98 /zm, separated by gaps of width 2 /mi, with period L = 100 /xm. 
To approximate the behavior of the resulting striated conductor, we apply the above 
X-array results. Using ([26]) with K s = 2H S , A c = d/2, 9 = nw/L = 0.497T, and 
tan# = 31.8, we find that the zero-field surface-barrier critical current of one of the 
narrow strips is I s o = 2H c i(%dLt&ii9) 1 / 2 = 2.9 A. The corresponding bulk-pinning 
critical current is I p = 2wK p = 0.98 A, and the bulk-pinning parameter fl32l) is p = 
0.22. Since p < 2/tt, the reduced critical current is i c = 1, and the self-field critical 
current of one of the strips is I c = I s q = 2.9 A, determined by the geometrical barrier 
alone. Since there are 50 such strips, the estimated total self-field critical current is 143 
A. The corresponding engineering critical current density, taking the cross section to be 
5 mm x 1 fim (ignoring the cross section of the substrate) is J e = 2.9 x 10 10 A/m 2 = 
2.9 MA/cm 2 , as opposed to the 1 MA/am 2 critical current density for the unstriated 
strip. The reduced critical current i c vs h for this case is similar to that shown in figure 
[7]for p = 0.2. However, the scaling field for i c vs h from (1781) is small: H sca i e = 5.7 Oe = 
450 A/m. To summarize, the X-array results using the above assumptions predict that 
the self-field critical current for the striated conductor can be significantly enhanced 
above that due to bulk-pinning alone (by approximately a factor of three using the 
above numbers). However, the application of a relatively small perpendicular magnetic 
field (of the order of tens of Oe or hundreds of A/m using the above numbers) can be 
expected to produce a strong reduction of the enhancement and to return the critical 
current nearly to that due to bulk-pinning alone. 
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Throughout this paper we have assumed that the bulk-pinning critical current 
density J p is field-independent. We now justify this assumption as follows. We see from 
figure [7] that the critical current has a significant dependence upon h = H a /H sca i e = 
B a /B sca i e only when h ~ 1. In the above two paragraphs we found for state-of-the-art 
YBCO coated conductors for which J p ~ 1 MA/cm 2 (10 10 A/m 2 ) that H sca i e < 10 Oe 
= 796 A/m or B sca[e = fi H sca[e < 10 G = 1 mT. Experimentally it has been found, 
for example in [54"1 [55] , that J p at 77 K in strong-pinning YBCO films is very nearly 
independent of field for H a < 100 Oe « 8000 A/m or B a < 100 G = 10 mT. Therefore 
it is an excellent approximation to assume that J p is field-independent over the range of 
applied fields for which the effects of geometrical barriers or edge pinning are relevant. 

In [53J, we argued that in state-of-the art unstriated YBCO coated conductors 
the pinning is due almost entirely to bulk pinning and that geometrical barriers (or 
edge pinning) have a negligible effect upon the critical current. The above calculations 
indicate that for striated coated conductors the additional edges produced by the 
striations could produce significant enhancements of the critical current in self-field, 
but these enhancements are strongly suppressed by relatively small applied magnetic 
fields. In applications where the magnetic flux density is of the order of 0.1 T or higher, 
where J p has a strong field dependence, we expect that geometrical barriers will have 
no significant effect upon YBCO coated conductors, even if they are striated. 
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